Nonlinear magnetoelastic behavior of the metastable bcc phases Co and Ni: 
Importance of third-order contributions for bcc Ni 
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The first- and second-order magnetoelastic coefficients of the metastable bcc phases Co and Ni are 
calculated by using a combination of the phenomenological theory of nonlinear magnetoelasticity 
with the ab-initio density functional electron theory. The magnetoelastic behavior of the bcc phases 
is drastically different from that of the corresponding fee phases. The recently synthesized bcc 
phase of Ni appears to be an example of a material for which third-order magnetoelastic effects are 
essential. 
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In recent years it became possible to stabilize 
metastable phases of materials by growth on appropriate 
substrates. For the transition metals Fe, Co and Ni this 
is especially interesting because in these systems mag- 
netism and structure are closely related. Using molecu- 
lar beam epitaxy, the fee phases of Fe and Co could be 
stabilized on substrates at room temperaturei. It be- 
came even possible to synthesize the bcc phase of Co£ 
on various substrates (Ref. ; 3 and references therein), 
and most recently^ the bcc phase of Ni. Both of these 
materials turned out to be ferromagnetic at room tem- 
perature, with a magnetic moment per atom of 1.53 /j.b 
(Co) and 0.54 /x B (Ni). 

From the viewpoint of technological applications of 
ultrathin magnetic films the most important feature is 
the magnetic anisotropy. Because in general there will 
be a lattice mismatch between the substrate and the 
magnetic film, magnetoelastic contributions to the mag- 
netic anisotropy may be important. For instance, it has 
been suggested^ that the in-plane anisotropy of bcc Co 
on GaAs is dominated by the magnetoelastic contribu- 
tion, although the epitaxial strains in this material are 
rather small, about 0.25%. For comparison, for bcc Co 
on Pt(001) the epitaxial strains are considerably larger— 
(-1.8% in plane and 5.1% out of plane). It is well known 
that for considerable epitaxial strains nonlinear contri- 
butions to the magnetoelastic energy become essential. 
This has been demonstrated experimentall y by cantilever 
bending beam experiments (see, e.g., Refs. I6l7ll8i) : When 
changing the direction of the magnetization in the epi- 
taxial film by changing the direction of the external mag- 
netic field, the magnetostrictive stress Act™ along the 
cantilever axis changes, resulting in a detectable change 
of the bending of the film-substrate composite. In the 
framework of linear magnetoelastic theory, this change 
should be independent of the magnitude of the epitax- 
ial strain and should be determined by the first-order 
magnetoelastic coefficients, i.e., B\ and B% for cubic ma- 
terials. Experimentally, however, a linear dependence of 
Acr m on the strain was found, which was ascribed to non- 
linear magnetoelastic effects. 

For a proof of this conjecture a knowledge of the first- 



and second-order magnetoelastic coefficients of the re- 
spective bulk material is required. The standard method 
to determine them is the ultrasonic pulse echo exper- 
iment. Because the attainable strains in these exper- 
iments are very small, it is, however, nearly impossi- 
ble to explore the second-order magnetoelastic coeffi- 
cients by these experiments. The first confirmation of 
the conjecture therefore was supplied by theory. By a 
combination of the phenomenological theory of nonlin- 
ear magnetoelasticity 9 with the ab-initio density func- 
tional theory it has been shown (see, for example, Refs. 
ll^llTlllT2lll^ and references therein) that the second-order 
magnetoelastic contribution indeed may be very large, 
especially for the case of Fe. The theory was also ableii 
to suggest a complete set of six cantilever experiments 
to determine the first-order (Bi and B-i) and the second- 



order (mj' z , m^^, m\ r , rn^) magnetoelastic 

coefficients of a cubic material. Thereby (mj' 2 , m^ 2 ) is 
related to pure tensile strains, (jn{ 2 , m^ 2 ) to tensile and 
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shear strains, and (ml' 2 , m^' z ) to pure shear strains. The 



first- and second-order coefficients have been calculated 12 
by the ab-initio electron theory for Fe, fee Co, Ni, NisFe 
and CoFe. 

The determination of the magnetoelastic coefficients is 
especially difficult for metastable phases which can be 
synthesized only as epitaxial films on substrates, like fee 
Co, bcc Co and bcc Ni, and in these cases the help of 
electron theory is very important. For the case of fee Co, 
the theory has showni^ that the nonlinear magnetoelastic 
coupling coefficients are essential for the magnetostrictive 
strain but have only little influence on the strain-induced 
out-of-plane anisotropy. In the present paper we apply 
the theory to the case of bcc Co and bcc Ni. It will be 
shown that in these systems the nonlinear magnetoelastic 
coefficients are again very large. Furthermore, it will be 
shown that bcc Ni is the first example of a system for 
which third-order magnetoelastic effects become relevant. 

According to Ref. ^] the magnetoelastic coefficients 
may be obtained by exposing the cubic material to cer- 
tain strain modes e,. Then the difference Ae^ in the total 
energy per atom when changing the direction of the mag- 
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netization from a\ to o3 is calculated: 

i = 1 : ci = e M = e , aj = (100), a 2 = (001) 
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t = 2 : e 2 = e yy = e zz = e , a] = (100), aj = (001) 



(1) 



Ae 2 = -Bxe + - (-Si - m 7 ' 2 + m^' 2 J eg (2) 
i = 3 : e 3 = = = e , a? = (010), a? = (110) 
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Ae 4 = S 2 e + -m i „ 
= e , aj = (110), of = (001) 
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The coefficient i?! and the pair (to 7,2 , to^' 2 ) of second- 
order coefficients are obtained from eqs.(l,2) by fitting 
parabola to the data points for Aei(eo) and Ae2(eo). 
Similarly, the coefficients S2, (m]' 2 , m^ 2 ) are obtained 
from eqs.(5,6) by parabolic fits. Finally, the pair (m^ 2 , 
TOj 2 ) is obtained from eqs.(3,4) via parabolic fits using 
the already determined coefficient to 7 ' 2 . As long as the 
parabolic fits represent the calculated data points Aei(eo) 
well, we can conclude that third-order magnetoelastic ef- 
fects can be neglected for the considered range of eo. 

The calculations of Aei(eo) were performed by ap- 
plying the ab-initio density functional theory taking 
into account the spin-orbit coupling which is responsi- 
ble for magnetoelasticity in a perturbative manner us- 
ing the second- variational method^ 4 . Furthermore, we 
use the WIEN97 codei^ which adopts the full-potential 
linearized-augmented-plane-wave method (FLAPW) 1 ^ as 
well as the local-spin-density approximation (LSDA) 17 
and the generalized-gradient approximation (GGA)i& for 
the exchange-correlation functional. The strains ti were 
applied with respect to the theoretically determined equi- 
librium lattice parameters a = 0.273 (0.281) nm for bcc 
Co and a = 0.273 (0.279) nm for bcc Ni in LSDA (GGA). 
The resulting LSDA (GGA) magnetic moments per atom 
of 1.63(1.74) Mb for bcc Co and of 0.47(0.53)/^ for bcc Ni 



are in agreement with the experimental values of 1.53 /ib 
and 0.54 ^b, respectively. 

For the case of bcc Co all the data points Aei(eo) could 
be perfectly fitted by parabola in the range —0.03 < eo < 
0.03, i.e., third-order effects can be neglected. Like for 
other materialsi2iii*i2iH the discrepancy between LSDA 
and GGA may be quite large. Because for the experi- 
mentally well investigated B\ of bcc Co, fee Ni and fee 
Co the agreement with the GGA values was better than 
the agreement with LSDA, we concentrate in the follow- 
ing on the GGA results. For bcc Co the values of B\ 
and B 2 are quite large as compared to bcc Fe, fee Ni and 
fee Co. The second-order coefficients are also large. It is 
interesting that there is a very large difference between 
bcc Co and fee Co. This holds even for the first-order 
coefficients B\ and B 2 which magnitudes are consider- 
ably larger and of opposite sign for bcc Co as compared 
to fee Co. In Ref. it has been assumed that for bcc Co 
the first-order magnetoelastic coefficients can be approx- 
imated by those of fee Co, in contrast to the results of 
our calculation. 

The case of bcc Ni is even much more interesting be- 
cause, as shown in Fig. 1, the data points for Aej(eo) 
show a drastic deviation from a parabolic behavior in 
the range —0.03 < eo < 0.03. To the best of our knowl- 
edge, bcc Ni therefore represents the first known material 
for which third-order magnetoelastic effects become very 
important. Another surprising result is that for bcc Ni 
the magnitude of B\ is very small (—1.3 MJ/m 3 ), much 
smaller than the one for fee Ni (10.2 MJ/m 3 ). As in 
the case of Co, the magnetoelastic properties of the bcc 
phase are drastically different from those of the fee phase. 
This is in line with the experimental observations^ that 
the cubic magnetic anisotropy constant K\ of bcc Ni is 
drastically different from the one of fee Ni, and this was 
attributed to the different electronic band structures as 
found by angle-resolved photoemission. 

We hope that our prediction of strong third-order con- 
tributions to the magnetoelastic properties of bcc Ni 
will initiate an experimental investigation by cantilever 
bending-beam experiments. To do this one has to grow 
epitaxial films of bcc Ni with various average epitax- 
ial strains eo which may be controlled with the film 
thickness^ and then the change A<r m of the magnetostric- 
tive stress due to a change of the magnetization direction 
has to be measured. For the case that third-order effects 
are relevant we expect a parabolic dependence: 



Acr m = a + £>ie + £>2e 



(7) 



As discussed above, a linear dependence has been ob- 
served experimentally already for several materials. The 
observation of a quadratic contribution would mean that 
for the first time a material was found for which the third- 
order magnetoelastic contribution is relevant. 
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Cl2 


C44 


A100 


Am 


fee Co±£ 


LSDA 


-15.9 


3.0 


243 


-53 


81 


102 


759 


796 


3.85 


2.26 




6.7 






GGA 


-9.8 


4.5 


184 


3 


59 


-41 


862 


1681 


3.13 


1.80 




4.9 




bee Co 


LSDA 


61.5 


-35.2 


-672 


575 


357 


-108 


-363 


336 


2.1 


2.95 


1.78 


48.2 


6.6 




GGA 


28.6 


-39.6 


-1013 


973 


148 


51 


-826 


611 


1.74 


2.33 


1.40 


32.3 


9.4 


fee Ni±£ 


LSDA 


12.6 


16.9 


-117 


23 


168 


-47 


-2 


388 


3.63 


2.20 




-5.9 






GGA 


10.2 


11.1 


-95 


71 


90 


-4 


108 


96 


2.95 


1.75 




-5.7 




bee Ni 


LSDA 


-0.2 


19.1 


-116 


-100 


917 


634 


« 


-308 


1.99 


2.62 


1.90 


-0.2 


-3.4 




GGA 


-1.3 


19.7 


-28 


-3 


87 


590 


-75 


-1347 


1.52 


2.32 


1.63 


-1.1 


-4.0 



TABLE I: The calculated magnetoelastic coefficients (in MJ/m 3 ), elastic constants Cu, C12 and C44 (in 10 11 N/m 2 ), the 
magnetostrictive coefficients A100 = — 2£>i/3(Cii — C12) and Am = — B2/3C44 (in units 10 ) Our calculated elastic constants 
for bec Co and Ni agree nicely with those given in Ref. Hsl 
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